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1.  Introduction 


In  this  report  we  present  a  series  of  medium-range  (10-dav) 
forecasts  performed  with  the  UCLA  general  circulation  model  (GCM). 

For  most  of  the  forecasts  we  have  used  more  than  one  version  of  the  GCM 
with  different  horizontal  resolutions  and  vertical  extents.  The  impact 
of  these  model  characteristics  on  the  accuracy  of  the  forecasts, 
therefore,  can  be  assessed. 

Section  2  contains  an  outline  of  the  UCLA  GCM.  Descriptions  of 
the  dataset  and  the  initialization  procedure  are  given  in  Section  3. 

An  inventory  of  cases  and  models  used  for  each  case  is  in  Section  4. 

The  major  part  of  this  report  consists  of  the  forecast  results, 
presented  as  maps  of  sea-level  pressure  (SLP)  and  geopotential  <Z) 
fields  at  selected  pressure  levels.  < 

2.  Outline  of  the  UCLA  GCM 

The  prognostic  variables  of  the  UCLA  GCM  are  the  horizontal 
velocity,  the  potential  temperature,  the  water  vapor  mixing  ratio,  the 
ozone  mixing  ratio,  the  surface  pressure,  the  planetary  boundary  layer 
<PBL>  depth  and  the  ground  temperature. 

One  of  the  unique  aspects  of  the  UCLA  GCM  is  the  treatment  of  the 
PBL .  This  is  assumed  to  be  well-mixed  and  is  made  the  variable-depth 
bottom  layer  of  the  model.  In  parameterizing  PBL  processes,  we  use 
bulk  assumptions  for  the  description  of  the  turbulence  (Suarez  et  al . . 
1983).  When  the  top  of  the  PBL  is  above  the  condensation  level,  it  has 
a  sublayer  of  stratocumulus  clouds.  The  parameterization  of  cumulus 
convection  and  its  interaction  with  the  PBL  is  based  on  Arakawa  and 
Schubert  (1974).  Details  of  this  parameterization  are  given  by  Lord 


(1978)  and  Lord  et  al .  (1982).  Surface  fluxes  of  sensible  heat, 
moisture  and  momentum  are  modelled  with  the  bulk  parameterization 
proposed  by  Deardor ff  <1972>.  The  entrainment  of  mass  at  the  PBL  top 
is  calculated  using  Randall's  (1979)  closures.  In  addition  to  the 
surface  drag,  the  model  includes  the  effects  of  vertical  redistribution 
of  momentum  by  cumulus  clouds  and  lateral  redistribution  by  a  non¬ 
linear  "eddy  viscosity". 

In  the  vertical,  the  model  is  based  on  a  modified  o-coordinate 
system  for  which  the  lower  boundary,  the  PBL  top,  and  isobaric 
surfaces  above  a  prescribed  pressure  level  (100  mb)  are  coordinate 
surfaces.  The  top  of  the  model  atmosphere  is  assumed  to  be  a  material 
surface.  The  vertical  finite-differencing  used  above  the  PBL  is  that 
described  by  Arakawa  and  Suarez  (1983).  It  guarantees  conservations 
of  the  global  mass  integrals  of  potential  temperature  and  total 
potential  plus  kinetic  energy  under  frictionless  adiabatic  processes. 

The  equations  are  horizontally  discretized  using  a  staggered 
latitude-longitude  grid.  The  scheme  for  the  horizontal  advection 
term  in  the  momentum  equation  is  based  on  the  potential  enstrophy 
conserving  scheme  presented  by  Arakawa  and  Lamb  (1981),  modified  to 
give  fourth-order  accuracy  for  the  ac vection  of  potential  vorticity. 

The  horizontal  advection  scheme  used  for  the  potential  temperature  is 
also  fourth-order  and  conserves  the  global  mass  integral  of  its  square. 
The  scheme  for  the  horizontal  advection  of  water  vapor  and  ozone  is 
based  on  the  same  formulation,  although  exact  fourth-order  accuracy  is 
abandoned  in  favor  of  a  formulation  that  prevents  the  occurrence  of 
negative  values.  In  all  other  terms,  including  the  continuity 


equation,  the  pressure  gradient  force  and  the  definition  of  absolute 
vorticity,  the  differencing  Is  of  second-order  accuracy. 

The  model  includes  both  diurnal  and  seasonal  changes.  Solar  and 
infrared  radiative  heating  rates  are  computed  following  the 
formulation  of  Katayama  (1972)  as  described  by  Schlesinger  (1976). 

The  calculation  includes  the  effects  of  absorption  of  solar  radiation 
by  o2one  and  water  vapor,  and  absorption  and  emission  of  long-wave 
radiation  by  ozone,  water  vapor  and  CO^  at  levels  below  30  km.  Above 
30  km  infrared  cooling  is  modeled  using  Dickinson's  (1973)  approximation. 
The  cloudiness  used  in  the  radiation  calculation  is  predicted  by  the 
model  by  assuming  overcast  conditions  whenever  the  predicted  relative 
humidity  exceeds  saturation  at  levels  below  100  mb  or  when  cumulus  clouds 
penetrate  above  500  mb.  The  parameterization  of  the  ozone  photochemistry 
is  described  by  Schlesinger  (1976)  and  Schlesinger  and  Mintz  (1979). 

The  CO2  amount  is  prescribed. 

The  geographical  distributions  of  surface  albedo,  ground  wetness,  and 
sea  surface  temperature  are  interpolated  from  prescribed  monthly  means 
based  on  the  observed  climatology. 

For  the  forecasts  described  in  this  report,  we  used  the  15-layer 
and  9-layer  versions  of  the  model  with  low  and  high  horizontal 
resolutions  (5°  in  longitude  by  4°  in  latitude  and  3°  in  longitude  by 
2.4°  in  latitude,  respectively).  The  15-layer  version  has  the  top  at 
1  mb  and  six  layers  above  51.8  mb  (with  Interfaces  at  26.8,  13.9, 

7.20,  3.73  and  1.93  mb),  resulting  in  a  stratospheric  resolution  of 
about  5  km.  The  9-layer  version  has  the  top  at  51.8  mb  and  is 
identical  to  the  15-layer  version  below  that  level,  except  that  the 


ozone  mixing  ratio  ts  prescribed. 
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The  GCM  has  been  evaluated  in  terms  of  long-term  simulations  of 
monthly  mean  fields.  Suarez  et_  al.  (1983)  presented  a  July  simulation 
with  the  low-resolution  9-layer  version.  A  more  extensive  evaluation 
of  that  version  in  terms  of  a  three-year  simulation  has  been  performed 
by  Randall  et  al .  (198S).  The  high-resolution  15-layer  and  9-layer 
versions  have  also  been  evaluated  by  performing  simulations  over 
several  months. 

3.  Data  and  initialization  procedure 

The  dataset  consists  of  the  FGGE  Level  1 1 1 — b  analysis  by  the 
European  Centre  for  Medium  Range  Weather  Forecasts  (ECMWF)  below  10  mb, 
and  the  U.S.  National  Meteorological  Center  (NMC)  global  height  and 
temperature  fields  above  that  level.  Winds  for  the  levels  above  10  mb 
and  latitudes  higher  than  20°  were  obtained  by  using  the  geostrophic 
relations,  while  for  lower  latitudes  they  were  obtained  by  linear 
interpolation  between  20°N  and  20°S. 

For  the  primary  prognostic  variables,  the  initial  conditions  of  the 
forecasts  were  obtained  from  the  dataset.  The  surface  pressure  was 
obtained  from  the  observed  sea  level  pressure  and  geopotential  fields. 
The  horizontal  velocity  was  then  obtained  by  interpolation  of  the 
observed  values.  The  geopotentials  at  the  interfaces  of  the  model's 
layers  were  also  obtained  by  interpolation  from  the  observed  values. 
This  particular  interpolation  procedure  is  designed  such  that  the 
model's  pressure  gradient  force  is  a  proper  interpolation  of  that  on 
the  pressure  surfaces.  The  potential  temperature  was  then  obtained 
from  those  geopotentials  using  the  model's  hydrostatic  equation.  No 

further  initialization  procedures  were  7  1 ormed  for  these  variables. 
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To  obtain  initial  conditions  for  the  ozone  mixing  ratio  and  the  PBL 
depth  and  those  for  the  water  vapor  mixing  ratio  with  some  degree  of 
consistency  with  the  model,  we  performed  an  "update  cycle",  in  which 
the  model  was  integrated  for  several  (generally  five)  days  preceding 
each  forecast  period.  For  this  integration,  the  initial  PBL  depth  was 
50  ml)  everywhere,  the  water  vapor  mixing  ratio  was  obtained  from  the 
dataset  below  300  mb  and  was  set  equal  to  2.5  x  lO-^  above  that  level, 
and  the  initial  ozone  mixing  ratio  was  zonal ly  uniform  and  taken  from 
Schlesinger  and  Mlntz  (1979).  During  the  update  cycle,  the  horizontal 
velocity,  the  potential  temperature  and  the  surface  pressure  were 
updated  at  regular  time  intervals.  The  water  vapor  mixing  ratio  was 
updated  by  preserving  the  relative  humidity  predicted  by  the  model. 
Below  10  mb,  the  updating  Interval  was  6  hours  during  the  FGGE  Special 
Observing  Periods,  when  data  are  available  with  that  frequency,  and  12 
hours  outside  those  periods.  Above  10  mb,  the  interval  was  24  hours 
since  the  OTIC  data  are  available  only  at  1200  GMT. 


4.  Case  Studies 

Tables  1  and  2  show  the  forecasts  presented  in  this  report  and  the 
versions  of  the  UCLA  GCM  used  for  each  case,  respectively.  A  detailed 
analysis  of  some  forecasts  is  given  in  Mechoso  et  al .  (198S,  1986). 


Table  1 


Inventory  of  forecasts 


Forecast 


Starting  Date 


1 

0000  GMT 

12 

January 

1979 

2 

0000  GMT 

14 

January 

1979 

3 

0000  GMT 

16 

January 

1979 

4 

0000  GMT 

21 

January 

1979 

S 

0000  GMT 

13 

February 

1979 

6 

0000  GMT 

IS 

February 

1979 

7 

0000  GMT 

17 

February 

1979 

8 

1200  GMT 

11 

June 

1979 

9 

1200  GMT 

1 

July 

1979 

10 

1200  GMT 

3 

July 

1979 

11 

1200  GMT 

5 

July 

1979 
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Table  2 


Versions  of  the  UCLA  GCM  used  for  each  forecast 
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